Decoherence and Loss of Entanglement in Acoustic Black Holes 
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We studied the process of decoherence in acoustic black holes. We focused on the ion trap model 
proposed by Horstmann et al. (Phys. Rev. Lett. 104, 250403 (2010)) but the formalism is general to 
any experimental implementation. For that particular setup, we computed the decoherence time for 
the experimental parameters that they proposed. We found that a quantum to classical transition 
occurs during the measurement and we proposed improved parameters to avoid such a feature. We 
also studied the entanglement between the Hawking-pair phonons for an acoustic black hole while 
in contact with a reservoir, through the quantum correlations, showing that they remain strongly 
correlated for small enough times and temperatures. 
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One of the main results of quantum field theory in 
curved space-time is the Hawking effect, i.e. the particle 
creation process that gives rise to a thermal spectrum 
of radiation outgoing from a black hole This effect 
together with its entropy complete the interpretation of 
black holes as thermal objects. On one hand, it is be- 
lieved that the heart of a theory that unifies quantum 
mechanics and gravity lies in understanding the nature of 
this thcrmality. On the other hand, it is also important to 
collect experimental evidence in order to gain insight into 
this phenomenon, but this is practically impossible since 
black holes' temperatures are less than nK. To circum- 
vent this obstacle, W. G. Unruh proposed an analogue 
gravity hydrodynamical model where phonons propagate 
in a fluid with a subsonic and supersonic regime [2| . This 
system obeys the dynamics of a massless scalar field near 
a black hole and provides an experimental implementa- 
tion to study the Hawking effect. Following this analogy, 
there have been several proposals that involved BEC [3( , 
moving dielectrics 0] , waveguides [B| , slow light systems 
Q, among others. 

The current proposals do not provide conclusive evi- 
dence of the Hawking effect (see for example Q). Nev- 
ertheless, we believe that a particular one by Horstmann 
et al, Q, provides a promising setup. This system con- 
sists of a circular ring of trapped ions moving with an 
inhomogencous velocity profile emulating a black hole. 
The signature of this quantum radiation is the correla- 
tion between entangled phonons near the horizon, Q, 
which can be measured by coupling the ions' motional 
degrees of freedom to their internal state, (lfjj . 

Given that we are interested in the quantum nature 
of the effect, in this letter we study how decoherence af- 
fects the measurement, which has not been considered 
previously. Moreover, we work with the field-theoretical 
description in order to present a derivation applicable to 
any implementation of an acoustic black hole. Neverthe- 
less, given our interest in the ion trap we specialize our 
results to that particular setup. 



The Model. We can describe the field associated with 
phonons in 1 + f dimensions, ^(x,t), by means of an 
action written in the suggestive form 



dtdx\/- 



(i) 



where the effective metric is of a Painlcve-Gullstrand- 
Lemaitre type, given by 



ds 2 



(c 2 - v 2 )dt 2 + 2vdxdt - dx 2 



(2) 



and v(x,t) is the local velocity of freely- falling frames. 
In the ion ring implementation, the ions have the clas- 
sical trajectories 9®(t) and the quantum field describes 
small deviations from them, and it is defined by ^(x = 
6i(t),t) = y/pS0i(t), where p is the mass density. We 
take the velocity profile as 
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where /3 = (u max + v nlin )/2 and a = (w max - v min )/2. 
The minimum and maximum velocities are related by 
the condition that each ion has to make one revolution 
during the period r. It is important to notice that we 
use an approximate velocity profile. The actual one must 
be smoother, but ours is a good approximation for our 
calculations. Taking this into account, the parameters 7; 
and 9h do not exactly match those of 

Although the problem can be treated in a discrete fash- 
ion, we choose to stick to the field description since the 
action given in Eq.([l]) is common to every implementa- 
tion of acoustic black holes, not restricted to ion traps. 
Then our procedure can be used to study the decoherence 
in any acoustic black hole. 

In this letter we propose to work out the non- 
equilibrium features of the interaction between the ion 
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trap and its environment, since these are responsible for 
decoherence. As usual for this kind of tasks, we use the 
Schwinger-Keldysh or in-in formalism. Our system is 
described by the action given in Eq.(fT]). Following the 
quantum Brownian motion (QBM) paradigm ll| , the en- 
vironment is described by a continuous array of bosonic 
harmonic oscillators, distributed in each position of the 
circular trap, and represented by the degree of freedom 
q v (6,t) with action 



Ss[q v ] 



dvl{v) \ d 2 x[ql(6,t)-v 2 ql(6,t)]. (3) 



The function I(v) is the mass of each environmental os- 
cillator. The interaction between the system and the 
environment is given by the action 



7 / dv \ d 2 x Hf(x)q v (x), (4) 



so that the total action is given by S'f^E', g^] = Ss^] + 
Ss Vh] +Sint Q" ] ■ F° r example, in the ion trap proposed 
in [8[ , the velocity profile is produced by an external elec- 
tric field, which is generated with surface-electrodes. Ir- 
regularities over their surface produce fluctuations and 
the heat bath models each mode of the fluctuating field. 
Details of the possible nature of the environment is ex- 
tensively discussed in [l2[ . The fact that we are thinking 
of a fluctuating field coupled with the coordinate of the 
ions justifies the bilinear coupling used here. To simplify 
the treatment we also assume an ohmic bath, although 
more general environment would not change the results 
13[. The bath is at rest with respect to the laboratory. 
Moreover, the reservoir naturally has the same discretiza- 
tion (Planckian) length-scale as the ions in the ring, and 
therefore it is also sub-Planckian. Certain effects devel- 
oped in the presence of a trans-Planckian environment, 
such as Miles-type instability 14J, are absent here. 
The elements of the density matrix are given by 
&~,q~\t) = (^ + ,q + \p(t)\y~ ,q~). We assume 



p(V+,q+; 



an initial uncorrelated thermal state between system- 
environment, both at temperature T = 7i(/cb/3) _1 . More 
general initial states would not change substantially the 
process of decoherence, see [l5|.The reduced density ma- 
trix, p I (^ ± \t), is defined as usual, performing a trace over 
the environment. 

We studied the evolution of the system following stan- 
dard procedure, such as [l^ . The evolution of the re- 
duced density matrix is derived from the effective ac- 
tion S cff [*+,*-] = - + Si F [*+,*-], 
where the index ± designates both branches of the time 



and the Feynman- Vernon influence action 



path, see 
is given by 



S w = J d 2 xd 2 x'A(x)(r>(x,x')Y 1 (x') + -N(x,x')A(x') 
where A(x) = - *~(a;) and S = 



^ (x)). Both the noise and dissipation kernels are lo- 
cal in space. Their exact expressions are the QBM ones, 
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D(t,t') = di>—!—smv(t-t')<d(t-t'), (5) 

Jo i\y)v 
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N(t,t')= / <Jv-j— coth^-cosi/(i-i')> (6) 
Jo i\ v ) v 2 

where ^ 2 jvl(y) = j 2 v for an ohmic bath. From S e s, 
we can write the semiclassical Langcvin equation of the 
field, given by 

-^=9^9^^)+ J dsT>(t,s)*{8,x)=t(x,t). (7) 

The field £ is a stochastic force with a gaussian probabil- 
ity density with zero mean and correlation (£(x)£(x')) = 
KN(x,x'). 

Decoherence Time. The problem with this approach is 
that the velocity profile depends on time in an awkward 
fashion during the collapse, making the master equation 
hard to derive. To avoid this drawback we first stud- 
ied the imaginary part of the effective action Im(S)[\E , c i], 
verifying numerically that the small period, named r c , 
when the velocity changes with time does not contribute 
to the decoherence time. We also show that the small 
angular region where the velocity is not constant is also 
irrelevant. Using these results, we will derive the master 
equation using the weak coupling approximation. 

For piecewise space and time independent velocity pro- 
files, the modes can be found using the usual null coor- 
dinates 

(8) 



t 



c(x) + v(x) 
d.v 



= t 



v = t+ I — =t-x v . (9) 
J c\x) - v{x) 

Each solution can be decomposed in u and v modes 
which have the form fy u i v (x, t) = ^> cos (wu/v). For 
fixed w, we studied the decoherence of each mode u and 
v, but we found that the result does not depend on 
this choice. Since the spatial dimension is compact we 
find the available uj demanding that \E r u / v (t, x + L) = 
ty u / v (t, x). This gives a numerable set of frequencies, 
and the maximum frequency is found by estimating that 
the wavelength cannot be smaller than the ions' sepa- 
ration, which we call 5 = L/N, with L the circumfer- 
ence of the trap. The steps towards the derivation of the 
master e quat ion are standard procedure, which can be 
found in 11, lldl. The diffusion coefficient of the mas- 



ter equation is d{t) = J Q ds cos(ojs)N(t,t — s) and we 
call V = J Q L dxcos 2 (ujx u / v ). The solution to the master 

equation is approximately p r ex exp{— T f Q dtd(t)} where 
T = 7 2 U(*+ - *-) 2 /2. The condition to find the deco- 
herence time is V dt d{t) rj 1. 

Since the proposal consists of measuring the entan- 
glement between phonons next to the event horizon, we 
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study the decoherence of two modes with the best reso- 
lution we could have to obtain an upper bound for to, 
taking (*+ - *~) ps y/pj. At T = 0, the diffusion co- 
efficient is, for an ohmic bath, d(t) ss wSi(wi) We 
also assume that we are able to see several oscillations 
of the mode before it decoheres, i.e. t ^> oj^ 1 . Then 
d(t) ps wtt/2 and io(T = 0) ps 2h 1 - 2 S' 2 {lottV)- 1 . 

In order to compute the decoherence time, we need an 
estimation for 7. In [§] , an upper bound for the stochastic 
force £ in units of the mean force on the ions F, £ = £F, 
was established, given by £ < £ max = 5 x 10 -6 , see Eqs. 
(80) and (81) of the second paper in FJ|. We relate £ with 
7 through the r.m.s. value of the stochastic force present 
in Eq. ([7]), using the £'s two point function. The relation 
between these variables is 7 = (y/2p/h(v max — i> m i n ). 
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FIG. 1. Plot of the decoherence time as a function of £, for 
maximum and minimum frequency, the velocity (right inset) 
and the system's initial temperature (left inset). 

FigfT] shows how the decoherence time depends on £. 
First, we check that uj max to ~ 10 5 3> 1 and the approx- 
imations used are consistent with the results. We note 
that for the stochastic force bound proposed in Ref.j§], 
decoherence is too strong, since the measurement lasts 
no more than r seconds due to classical instability issues 
of the system, so to must be bigger than t in order to be 
able to measure quantum Hawking radiation. Let's say 
that the decoherence time must be a couple of orders of 
magnitude larger, to ~ 100r, then £ must be less than 
£max = 2x 10 -8 . We also show the dependence of the 
decoherence time on the velocity profile in FigfT] We plot 
to as a function of u m i n , for £ ma x, concluding that it does 
not change substantially since the dependence is smooth. 
Of course, for low enough velocities to can be less than 
t, since it is monotone. Nevertheless, w m i n must be close 
to (27t/t)0.83 in order to produce an event horizon. 

For the case of small temperatures, we separate the 
v integral between two intervals (0,2/3 _1 ) such that for 
small enough frequencies fiv/2 <C 1, coth(/3i>/2) ps 2/fiv, 
and for high enough frequencies (2/3 _1 ,oo) such that 
^> 1, one can take coth(/3z;/2) ss 1. For times 



longer than the period of the mode, wi > 1, the dif- 
fusion coefficient has the lower bound J Q D dt d(t,(5) = 
cunt ol '2 + 4w~ 2 /3~ 2 , resulting in a decoherence time of 



h 2 



Vmin) 8 2 U1Tp 2 V uj 3 tt[3 2 



o(/r 3 )- 

(10) 

In the FigUJ it is shown the dependence of the decoher- 
ence time with temperature for several couplings. For 
£ ~ 10~ 8 , the decoherence time is much larger than r for 
a wide range of temperatures. Therefore, it is clear that 
our proposed coupling is a significant improvement from 
the one proposed in |8j , and composes the most adequate 
parameters for the desired experimental conditions. In 
[8| it is proven that it takes some time to generate the 
entanglement between the Hawking phonons, even in an 
equilibrium situation, and that the entanglement is only 
appreciable for temperatures below IOOTh- Taking this 
into account, and even though it takes a large temper- 
ature to affect the quantum features of the system with 
such a small coupling as £ max , our results are only useful 
to measure the Hawking effect below lOOTff . Moreover, 
it can be verified that the low temperature limit is only 
valid below - lOOTjf. 

Correlations. Hawking radiation can be understood 
as a pair production of virtual particles, one of which 
falls into the black hole and the other, outgoing, be- 
comes real, building up the Hawking radiation [l|. The 
role of the correlations between this pair of particles has 
been studied as a signature of the quantum Hawking ef- 
fect, for example It was found that the entan- 
glement between this pair is translated to a sharp peak 
of (II - (xi, t)H~ (x, t)) as a function of x, where n is the 
canonical momentum conjugate to VE' and n~ corresponds 
to its left-moving modes, x\ is inside the black hole and 
x is outside. This magnitude was calculated using the 
Israel- Hartle-Hawking state, and it was also calcu- 
lated numerically in the case of the circular ion trap, 0]. 

We derived this quantum correlation and obtained an 
analytic expression that depends on time and also on 
the initial temperature of the system. We also show 
that before decoherence the difference between the closed 
and open system correlation is indistinguishable, but for 
later times the correlations change substantially from the 
closed case one. To obtain the correlation we follow the 
procedure presented in [l9[, Eq.(4.8), which consists of 
solving the Langevin equation, multiplying the solutions 
to construct the correlation and then integrating over the 
initial state and the stochastic source. Since the system 
is weakly coupled to its environment we treat the dissipa- 
tion kernel term perturbatively, calling ^ and G ro t the 
solution of the free wave equation and its Green function; 
thus the full solution is * + G Tet ■ (£ - D • 

We first obtain the analytic expression without en- 
vironment. For this we first compute the left-moving 
modes in a background, given by Eq.(f2]), which satisfies 
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FIG. 2. Plot of the momentum-momentum correlation for 
Xi — — 10a and t = 100r c . The profile parameters are 
fmin = 0.9, Wmax = 1.1. It is shown for different tempera- 
tures. Inset: Relative contribution of the correlation due to 
the environment as a function of time, in units of the collapse 
time t c . 



due to the environment becomes of the same order as 
the correlation of the closed system. We expect that af- 
ter this happens, the correlation gets distorted, the peak 
gets lost and the entanglement between the Hawking pair 
gets lost due to decoherence. 

If an experimental proposal is to be taken seriously, 
then decoherence must be taken into account. We hope 
that this study reassures that the measurement of the 
Hawking effect is possible. We also provide a deriva- 
tion of the correlation between the Hawking phonons 
as a function of time and temperature and we check 
that the relative contribution of the environment to this 
magnitude is irrelevant. We hope that this would help 
to deepen the understanding of the Hawking effect in 
acoustic black holes, and the entanglement between the 
Hawking-pair. 
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8l^~ = 0, where 8l = (dt + vd x — d x ). The profile 
is v(\x\ < a,t) — cr(t)(l — kx), v(x > a,t) = <r(t)v m i n , 
v(x < a,t) = cr(t)v max , where er(t) is a function that sat- 
isfies ct(0) =0 and a(t S> r c ) = 1, and t c is the time of 
collapse of the acoustic black hole. We solve the equa- 
tion using the method of characteristics, well known from 
fluid dynamics of compressible fluids [20( . We normalized 
the modes in a way that resembles the usual Minkowski 
expansion at t = 0. 

When solving by the method of characteristics we get 
two solutions; for example, we get the curves in the x > a 
region and the x < a region. The curves originated in 
xq < a must match along the boundary x = a. These so- 
lutions present the Hawking radiation-related peak. Nev- 
ertheless, when the characteristic starts from x — a or 
bigger, then the full solution is the v(x > a,t) = cr(i)w m in 
one. In this case the peak is not present. Therefore, the 
characteristic with xq = a provides a boundary between 
the region where the horizon is important and the region 
where it is not. For the regime where the Hawking cor- 
relation is present, the result is plotted for different tem- 
peratures in Fig 121 As noted in the correlations are 
present even when the ions' temperature is higher than 
the Hawking temperature, and this seems to be a general 
feature of the Hawking effect, as long as decoherence is 
not present. 

We also estimate the relative influence due to the pres- 
ence of the environment, e r = \(C c (k) — C (k))/C c (k)\, 
where C c is a single mode of the correlation (Hi^) when 
the system is closed and C is the correlation when the 
environment is turned on. In the inset of Fig[2]wc plotted 
e r as a function of time, for low wavenumbers, e r (k ~ 0, t) 
since it decreases with increasing k. For small times, 
e r ~ and then the environment does not affect the cor- 
relation. For large enough times, the relative correction 
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